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ABSTRACT 


The Sei ico vv resistance and self-propulsion model test results are 
Mare senteda by curve fittins methods wnick y==nit in functions containing 
Chebyshev polynomials. The accuracy of the curve fitting methoa is dis- 
Eu, and the range of validity of the functions is presented. 


EZconpUu5er program for the calculation of EHP and SHP for. single-screw 
merchant Ships with null proportions falling within the range of the 
Series €0 is discussed in detail. Several examples in the use of the pro- 
gram are presented. 

۱ ۶ concluded that estimates of EHP and SHP can be calculated to 


within a tolerance in the neighborhood of 1% on the average, and with great 
certainty, to within a tolerance of less than 3% for any given speed. 
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I. INTRODUCTION 


1. Foreword 

Tne determination of the required horsepower necessary to fulfill 
beer speca requirements of a new ship design is one of the many problems 
which faces the naval architect in the early stages of ship design. 

He must have some way of knowing what effect various choices of size, 
Meepertions, and fullness of the ship will have on the power required to 
propel the ship at the speeds he desires. | 

In determining this power requirement, the naval architect can draw 
upon his past experience in ship design and the results of specific isolated 
model tests which are published periodically in the various technical 
Journals of the profession. 

Ever we results of experiments conducted on families of models; 
Ica the different design parameters are varied systematically, are the 
naval E ce’ most valuable sources of reference in estimating the 
mower required for his new design. The result of these tests are usually 
eepressed in the form of design charts from which the naval m 
using interpolation techniques where necessary, can calculate the power 
requirements for a number of combinations of design parameters, and thereby 
mec an intelligent choice of the best combination of parameters to give 
imum power with due consideration for the other limitations of his design 
conditions. The Series 60 is such a family of models designed to aid the 
naval architect in the estimation of power requirements for single-screw 
merchant ships. 

En The Series 60 Model Tests 
The research on Series 60 was carried out at the David Taylor Model 


Basin (DTMB) of the United States Navy. Experiments were conducted as part 


E3 





of the Bureau of Ships Fundamental Hydromechanics Research Program during 
the years 1948 to 1960. The results were published in a series of papers 
presented before the Society of Naval Architects and Marine Engineers, 
and ۲ final report, incorporating all of the information contained in the 
interim reports was published in July of 1963, (1)* 

erence survey Of current practice in the shipbuilding industry to 
ensure, as far as possible, that the series would cover the normal range 
of proportions of modern ships, proportions were chosen for five parent 


models of the Series 60. These proportions are listed in Table I below. 


TABLE I 


Hull Proportions for Five Parent Models of Series 60 


| | | |> |> |= 



















۰ 
T 
i 


aft 1.54 aft 1.5 2.5 Ta 


s Orne block coefficient and LCB position, a total of nine models were 
tested in which the L/B and B/H ratios were varied, for a total of forty-five 


models. To cover the general spread of L/B, B/H, and A /(1/100)? for 





+ 
Numbered references throughout the text are listed on ra 


EO 





existing designs and the possible variation in LCB position, & grid was 


adopted as shown Sy the dottei lines in Figure I. 


x 
FIGURE I 


۳ ۰ 2 ` 5 
Variation of Proportions with C for Series 60 Model 


> 23 
~ ~ 1 










* NEW PARENT MODELS 
FOR $ AN0 Ñ VARIATICNS 





0.60 0.65 0.70 0.75 0.80 0.85 
BLOCK COEFFICIENT Cg 


u 





* From reference (1), rege II-3 





The pattern for a typical case (a, = .60) is shown in Figure II. 


x 
2 1 ٦۳ 


Typical Variation of L/B and B/H Ratios for 


a Given Value of G, 





6.5 7.5 8.5 
+ 
B 


s (pip Section of the models had no desdrise in accordance with 
Merenu practice. The relationship between Cr and C for the Series 60 is 
۱ P 


as follows: 


C, = 0.0857 Co + .925 
End, since Cp = Co x 
2 
Œ = 0.0857 C + ۰.925 C ii 
D 2T C, 225 C, (E) 


LWL for the Series 60 hull forms is related to L by the following 
simple formula: 


WL = ho6.T5 x L — |. (2) 
400. 


* From reference (1), page II-3 


he 





The results of the Series 60 resistance tests are presented in the 
form of easy to use design charts, and in numerical roint-by-point form 
LED Les. 

Values of total resistive coefficient تا‎ for forty-five ships of 
length between perpendiculars of LOO feet, and varying values of و‎ BB, 
and B/H are presented in Tables B-1 through B-45 of reference (1). 

Utilizing the Froude assumption that the total resistance can be divided 
into two parts, the skin friction of an "equivalent plank" and the 
residuary resistance, the values of C include a skin friction 6967 ٦ 
(C,.) calculated for salt water at 59°F in accordance with the American 
Towing Tank Conference (ATTC) 1947 eee and an additional allowance for 
ship correlation (AC 3 of + 0.0004 as recommended by the ATTC 1947. 
Methods for scaling the values of C. for an arbitrary ship lens Mm be 
Be eussed in Chapter IV. 

Values of propulsive coefficient (PC or EHP/ SH?) Tor forty- 086 
of length between perpendiculars of 600 feet, and varying values of وا‎ 
L/B, and B/E are also presented in Tables B-1 through B-N5 of reference (1). 
Unlike the values of ون‎ tne propulsion data can not be corrected for 
variation in ship length. The choice of a 600 foot ship to illustrate the 
propulsion tests was mace principally because it was considered more rep- 
resentative of modern ships than the 400 foot length chosen for the resistance 
presentation. * The propulsion deta presented was obtained using a standard 
propeller diameter equal to 0.70 of the designed load draft. 

Pata for the effect of variation of LSB position on power requirements 


is by no means complete. Corrections to total resistance for ٦37307 


Also known es the Schcenher line. 


xx 
Reference (1), page VII-k. 





LCB position for five values of ° are presented in Tables 49 through 
53 in reference (1), with only four or five different positions of LCB 
noted. 

Although the propulsive data is presented for ships of 600 foot 
length, if the designer adopts the lines of Series 60, the position of 
۳ Used in the parent forms, and a propeller having the standard 
ratio of diameter to draft of 0.7, he can make a very accurate estimate of 
ES tne EHF anû SHP of a ship for any particular selection of L, B, B, 
and A. 

3. Previous Work in Series 60 Computerization 

iiescevelOpment Of a computer program for estimating EHP using one 
Series 60 results was first ت5‎ 6 by D. Khousry in 1962. 9). "his 
attempt represenicu Liz “ries 60 by curve fitting methods to some extent 
meee btca surface, but it relied exclusively on interpolative methods 
Bor C. wnich are both cumbersome to program and highly consumptive of 
Eornputer storssze space“. Purthermore, the large amount of storage srace 
megurred in this tyre of program made it unsuitable for use as an EHP or 
SHP subroutine in a larger program. It also should be noted that an 
Per oolacive scheme is not very efficient since tae test Gesuits are 
mot presented in & uniform manner, i.e. at low values of Co models were 
tested up to a V/VLÄL of 1.10 whereas the higher values of e. were vestel 
up to a NEE of 0.80. This problem was avoided in the Khoushy work by 


discarding a large amount of Series 60 data in the low and high ranges of 


V/YLWL. 


x 





The storage of approximately 6000 numbers is required for the Series 60. 
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Therefore, it was decided that curve fitting would be the most 
advantageous method of representing the Series 60 model test results 
fee computer program. 
lh. The Curve Fitting Approach 

The curve fitting approach for nen NUS Series 60 results 
EN rv desirable since Once a function is found, which represents tne 
ven date, one knows the error he can expect in predicting EHP and SHP 
۰” 16 6461۷606 functions. The fitting of curves through the Series 
60 resistance data was started in early 1966 by T. M. Pitidis-Poutous. 
îê theory of the curve fitting method employed was developed in 
reference (3) in 1965 with the mathematical derivation of ships lines 
meine tne focus of attention. ‘Tnis curve fitting method, as it applies 
to the representstion of the Series 60 resistance data and its extension 
Eu esent the propulsion results for this work, is described in 


Chapter II. 





II. DESCRIPTION OF CURVE FITTING METHODS USED FCR THE SERIES 60 


The theory of the curve fitting methods used for representing the 
Series 60 model test results stated herein is largely adopted from reference 
(h). 

Since the chart reading methods of the Series 60 are not suitable for 
Bemputerization, it was decided to fit curves through the numerical results 
presented for A and PC (or EHP/SHP) in Tables 3-1 through B-45 of reference 
B. 

In the Series 60, both C, and PC are functions c four independent 
variables, i.e. V/YLWUL, رت‎ L/B, and B/H. Thus, either C, or PC can be 


represented by an equation of the form: 


Mz 


3 


R:l 


PE 


2 


4 


Oak i $ (5) y. (3) X (C) o. تا‎ (3) 


=> 


3 
"N 


i 


C. 


Semone problem reduces to determining the matrix of coefficients ] 0 il 
from the given data. If the coordinate functions $ ۰ Y , X, and NQ are 


ken as Ordinary polynomials, the above expression becomes 


I J 


Fr Y» O aki (sy [C 309 t رت‎ 


atl Jz! اک‎ nel 
It can be shown that there are certain advantages to be gained if the 
coordinate functions are taken as the Chebyshev polynomials, with no 


개 | : * 1 EA 
additional complications. The independent variables V/ PLL, D L/3, 


* See reference (5). 


۳ 





end B/H can be normalized in the range of -l to +l by the equations given 

MN Teble Ii, and the corresponding symbols used for the original, normalized, 
EN 56er prosram variables are given in Table III. The equations of 
Table II were derived on the basis of the limits of the Series 60 model 


tests. : 


Stipulations on the maximum powers I; ول‎ Ky and 5, to be used in 
equation (3) are determined by the amount of data available from the 
Series 60 results. For example, each model was tested for three different 
values of B/H and L/B, therefore I and J must be equal to or less than 3. 
Models were tested for five different AoE eG C therefore K must be 
equal to or less than 5. Thus, the number of terms in the equation will 
be equel to orless than 45N. The value of N cannot be arrived at in the 
Ame vay as 1, J, and A, Since each model was tested for twenty to thirty 
different values of VEL, یت‎ ۶۲67۰770۶ irom model to model: 77 
for the moment, that N = 10 would be adequate, we see that the determina- 
tion of the coefficients Any; ; would involve the solution of 450 
linear alsebraic equations. This is virtually impossible, both from the 
numerical and computer storage points of view. 

These considerations prohibited the use of a direct, 4-parameter fit, 
wnai cated tke necessity of a step-by-step technique. 

Using the symbols for the normalized variables as defined in Table IIT, 


equation (3) becomes: 


III a: EAN, m) 


oT 


y 


it 


A =| 


2 2 Ane Rr) 1, (7%) (h) 





TABLE II 


Equations for Norralizins Indevendent 
Bu! Te ne, INES DENT 이 


Karsanles in Series 60 Range 





E25 — - 1.875 N 77+ 
== 





ELLOS 210 ; L = 1.5u + 7.0 
1.5 B 

No P - 6.0 ; = = 0.52 + 3.0 
H H 

m= 10.0 C - 7.0 ; 0 = 0.lx + 0.7 
b b 


2102 





TA 


20 






10 


Hd 


III ۵ت‎ 


Wod 


WLO 


gTOoquÁg ue13014g 3۹700100 pDug ‘pəz F ] 01100 T BUTS TIO Jo 3011301018 90000 


Toquis 593000 


3011001 ]۲ 92 و۲00۳ 


0 و۲00۳ 


EI 





where: NOS JT(<-1) + IT CE 


and Fo = à 2 (a) Y (3) Dt (x) 


The above substitutions reduce the four-dimensional matrix "a" to 
a two-dimensional matrix "A", which is easier for computations, and 
femmes rhe problem more appropriately for the method to be used. For 


given values of u, 2, and x, equation (4) can be written as: 


N : 
n ۳ > 20 9» (+) oe 
with l L F... | 
<a z > 0 7 ( 4,5, X) (6) 
2:۱ 


Equation (5) represents the variation of Y with respect to v for 
the pth model and the N coefficients (C54; I 2 N) can be evaluated 
from the data in Tables B-l through B-45 of reference (l), using the 


least-squares criterion: 
2 


N 
a (Y): 7 2 Con Oo (vx) = MINIMUM 


4 
by differentiating partially with respect to the unknown coefficients. 


ode Ss pere titted in this fashion, and the resulting equation 
Tor each model was evaluated at intermediate values in order to ensure 
that the polynomial had no undesired oscillations. In fitting the curves 
FOr وت‎ it was found that a reasonably good fit was obtained with N = 10 
(corresponding to a 9th order polynomial) and since this is rather high, 
۳ 95 19 (1 ۳ Ol simunogucing a Suitable weighting function was explored. 


TOE MD ICID Ee repeated with the following two alternative 


BIO. 





Weighting functions assigned to oe 


Y (c, * re*) ( V/vzaz) 


Y (C2x 107%) / ( vis y 


The last alternative was found to give the best fit with N = 10 for most 
models. Going from N= 10 to N= 11, the maximum error experienced in 
C, x 10? only decreased from .025 to .024 for Table B-1 of reference (1). 


Having obtained matrix Cc.) corresponding to equation (5) with 


(number of models), matrix | was obtained by exact fit, i:e.‏ 45 = وو 
L= 45. It should be mentioned here that a certain amount of regression‏ 


pauta also be tried on this step of fitting, and this was actually 

E supred. but without success. The reason is that, although it might 
be possible to obtain a good fit hor = or PC at a constant value of 
V/TLUL with L less than 45, it does not follow ae such a fitting could 
be applied successfully to the coefficients of the C, - V/YLiL equation 


(i.e. C Thus, direct h-parameter fitting of the data could have 


pa) 
produced an equation with less error and less terms, if it were possible 
Bumerically. 

Tne values of PC were fitted in the same manner described above, but 


no weighting function was used. 


The final equations obtained are as follows: 
¡O ےھ‎ 


C, x 10° _ w A E x 2 (v 
(v/ Vis) 22 i un | s ۱ 





lo 4$ | | 
1 2 PE. D, Fe (4,9,*) TS (v] 
72.1 ۸ < ر‎ l 


2 





where 


Ie TB Ta ۳ 
J + 3 ) 2۱٦ + 2 ( k-1) 


Br = P t. Chebyshev Polynomial 


z 


>. 
1 


The elements of the matrices (A ¿Jana (3 n are read in as data 
n n 


ma tne main program, and are listed in Appendix A. These matrices forn a 
compact computational package of 900 numbers for the calculation of EN 
ISO ft. ship) and PC (600 ft. ship). 

No attempt was made at مت‎ curves through the cata for correction 
mo lots! resistance due to variation of LCB position. The inadequacy 


۳۳۲۲ ۱ date presented, as well as the erratic nature in which it behaves 


Beeren, »recluded the use of the curve fitting methods described above. 


Eh 





III. RESULTS OF CURVE FITTING METHODS USED AND RANGE OF VALIDITY 


The results of the curve fitting methods used for C. and PC were 
first tested by reproducing the data contained in Tables B-1 through B-45, 
i.e. for the forty-five models of the Series 60. The results of this 
i cure presented in Table IV, where the maximum error, in percent, 
ENDerienced in C. and PC is noted for each of the forty-five models. 

The maximum error experienced in ER was 2.740% for the parent model 
1218. The maximum error experienced in PC was 2.030% for model 4276. 
Although these Tires appear prohibitive at first glance, it should be 
emphasized that these are the maximum errors experienced for any 7٣ 
model, and they occur at only one value of V/VLWL. The errors at other 
values of آ۷‎ for any given model are well under 1% in the great 
때, or cases» In fact, as can be seen from Table IV, the maximum 
ors experienced in Ce and PC for any given model are well below 1% 
En most cases. To illustrate Chis point further, plots of e versus 
TL are presented in Figures Ill through VII, and plots of FC 
(or EHP/SHP) versus V/VLWL are presented in Figures VIII through XII, 
for the five parent models of the Series 60. In these plots, the solid 
line represents the value " C. and PC as calculated from equations (7) 
and (8) respectively, and the points represent the value at that V/YPLWL 
presented in the tables of reference (1), which are in error enough to 
be visually discernable on the scale used in making the plots. It is 
۱ 6 ۳۵ هه‎ Faaty Figure IV represents the worst case experienced in 
W ceing either ^ Or re. 

As a second test, equations (7) and (8) were evaluzted by generating 


O 


125 interpolate and extrapolate models, especially concentrating on 


ALA 





TABLE IV 


Maximum % Errors Experienced in C, and PC for the 45 


¿ana Fe tor tne 7) 
Models of the Series 60 


۴ Denotes parent model 


Series 60 Model Max. % Error Max % Error 
Table No. _No. C, IE BE _ int, — in PC 
B-1 Io «So 6.5 25 .578 .853 
B-2* OG 6O 7.5 2.5 1.230 . 618 
B-3 mer nes 8.25 2.5 1.515 . 672 
po) ho6h .65 6.25 2.5 . 0 .616 
B-5* ۰۰٠ 65 7.25 2.5 일 0 .298 
B-6 EO 2.5 , 701 . 696 
B-7 I0 6.0 2.5 .531 . قبام,‎ 
B-8* ISO 7-0 2.5 و‎ 402 
B-9 ۱ 0 8.0 2.5 2202 B 
۲۶-0 4268 .5 5.75 2.5 "O19 .386 
B-11* ۳۳ 75 6.75 2.5 poe 697 
B-12 1۳۰۰٦٠٠٦٠ 7.75 2.5 1.045 .682 
B-13 ION SO 5.5 2.5 01 .372 
B-14x Wee OO 6.5 2.5 .381 .134 
B-15 Be tas 2.5 243 .394 
B-16 ۱۵52 .60 6.5 3.0 .885 437 
B-17 B8 55 7.5 3.0 . 714 . 393 
B-10 ho5h .60 8.5 3.0 Su 1.261 
B-19 DOSES T 6.25 3.0 . 602 . 529 
B-20 ۱۳9 سا‎ 7.25 3.0 .1405 46 
B-21 HO) 8.25 3.0 .600 .516 
B-22 500. 0.0 3.0 وم[‎ TT 7 
B-23 ٠٦0 7.0 3.0 zu .590 
B-21 EDO JO 8.0 0 278 ae 
B-25 2۰۰۰۰۱۰۰۶۰ 6 0| 2528 2.030 
B-26 129 ۰.75۶ 6.75 3.0 ۳ 1.430 
B-2 T o ml T.T 3.0 . . 291 
B-23 ۱ OO 5.5 0 ‚oh, 07 
B-29 1509022950 6.5 3.0 16 .129 
B-30 ۱۱2۸62 ۰.0 7.5 3.0 .082 20h 
B-31 2 600 6.5 3.5 ~ 54) lio 
B-32 h253 .60 7.5 3.5 01 878 
B-33 은 8.5 3.5 59 .906 
B-34 ۱۳ ری‎ 6.25 3.5 6344 37 
B-35 DONNEES 7.25 3.5 u . 359 
B-36 4066 .65 8.25 3.5 . 493 . 661. 
B-37 ٣٠ 6.0 3.5 .203 , 3 
B-33 ۷0 3.5 .135 n» 
B-39 ۱۳ 3,5 312 , بابا‎ 
B-40 PONES US 3.5 ۰.223 .598 
B-41 DNE CCS TS 3.5 . 506 se 
B-42 O 3.5 .919 1.275 
B-43 ۰ ہہ‎ 1 5 .192 2 
B= H1 neo 80 166 3.5 .191 SS 
B-45 1250 .80 7.5 3.5 2 .303 
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intermediate values of the fiveCb's used for the Series 60 models. This 
Nas done in order to determine the range of validity of these equations 
Mem respect to their four independent variables. It should be emphasized 
that the ranges of validity determined below are not just a result of 

the inadequacies of the curve fitting methods used; they are very much & 
۰۰۰۰7 7ہ‎ the limits of the Series 60 model tests Eee and the 

range of data made available in reference (1). 

The data presented in reference (1) is based on the use of three values 
of B/H, i.e. 2.5, 3.0, and 3.5. Since equations (7) and (8) rely ona 
second order interpolation scheme in B/H, which is the method suggested 1 
reference (1) when using the design charts, there should be no discrepancies 
caused by this variable. It is es oud that € limited amour 
Abra poletion is possible in s/z." 

Similarly, there should be no discrepancies due to L/B since equations 
(7) and (8) rely on a second order interpolation scheme for this variable, 
which is also the method used by DTIB in the preparation of the Series 60 
Mesien charts. However, as can be readily seen from Figure I, there are 
limits (imposed by the proportions of the Series 60 models) on the range 
of interpolation possible in L/B. The range of interpolation in L/B isa 


function of C and can be expressed as follows: 
p 


MM + 9.0 - 10 Cr (9) 
2 

이 1 0 1 0 - 10 C, (10) 
2 





¥ 
See example of the ship SCHUYLER OTIS BLAND in Chapter V. 
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A small amount of extrapolation beyond eee limits should not cause 
broubie. 

In producing the Series 60 design charts, & third order interpolation 
in C was used by DTMB. In equations (7) and (8) a fourth order inter- 
Bolative scheme for E. is used,which could produce slight discrepancies 
En interpolation and extrapolation. It is recommended that the program 
be used only in the range of ° OCCO bo ٦ 

The Series 60 model tests impose definite limits on the range of 
V/VLWL permissible. For example, in testing equations (7) and (8) with 
the 125 models mentioned above, it was found that with C = 0-775, B/E = 3.0, 
and L/B = 7.0, the values of C, end PC given by equations (7) and (8) 
dropped off sharply at V/VLUL = .925 end esco mont negative at V/ VLL = 950. 
MS Ss understandabic Site at these values of V/TENL we ace outside the 
range of results given by the Series 60 tests. Therefore, to avoid 
Mecurrences of this nature, a ‘trap’ has been built into the EHP-SEP program 
Meceeri0ed in Chapter IV, which puts limits on cae va tuoo V/YCTL allowed , 


depending upon the value of 5 These limits are presented in Table V. 


TABLE V 


Limits of V/Yi;L Alloved by Program 


Pange of C 
nge of C, 








Within the ranges of independent variables stipulated above, one can 
expect values of C, and PC from equations (7) and (8) to agree with the 


published results of Series 60 to within 1%, with consistent reliability . 


아스 





IV. DESCRIPTION OF THE SHP-SHP COMPUTER PROGRAM FOR SERIES 60 


| The program for calculating EHP and SHP using the Series 60 model 
test results was written in che FCRTRAN IV language. A flow chart for 
ihe program, a listing of the FORTRAN statements, fama the 900 elements 
of the matrices کا‎ ١ 83 (3.,) mentioned in Chapter 11 are included in 
En sues A. The following is a general description of the methods of 
Me lculation used. Detailed instructions for the use of the progran are 
contained in Appendix C. 
ES Inputs 

Meese one elements of matrices (a) and ta are read nm + c 
The number of ships to be calculated (NJOBS) and an identifying alphameric 
EE 1s then read. It takes three cards to dessribe eäch ship to 
be calculated. The first card descrites the irvut mode (INOPT) for the 
EE DOLO be Calculated. The operator is afforded an option of one of two 
possible input modes for each ship. 

@. “Kode I 

Entry is made with L, B, H, C s Vmax» Vmin» Vinerement:» AC. Position 
of LCB from y y S, P, and kinematic viscosity. 

۲ rode II 

Entry is nade with L, B/H, 4, Cy, Vnax» Vmin» Vincrement» 4Cp, Position 
of LCB from % j رص رت‎ and Kineratic viscosity. 
2. Calculation of necdedzarameters 

Using the input information, the program then calculates the remaining 
parameters required for executing the remainder of the program. For example, 
if the eqs rede using Mode I, the program calculates 5, B/E, Co using 


equation (1), the Series Ó0 equivalent LWL using equation (2), the volume of. 


E 





displacement V , midship coefficient C,, and L/B. 


3. Normalizing hull proportion parameters 
The program calculates the normalized hull proportion parameters 


for L/B, B/H, and C, using the equations in Teble II. 


A. HD to make sure ۷/11 is within range of model test 


The program calculates V/YLWL maximum and minimum and ensures that 
| these values are within the range of validity given in Table V. If they 


are not, the program alters V i, and V .. accordingly, until they are 


n 


thin the range. 


C Iculation of number of increments in V to be calculated 

Using Virax> Vine @2° Vincreament» the number of different speeds (m) 
Bench EHP and SHP are desired is calculated. There is @ maximum of 
u (reni speeds allowed for each ship. If this number exceeds 100, 


"m is 100 O Mess Tf the 


the program keeps doubling Vincrement until 


Mera tor desires EHP-SHP information for only one speed, m is set at one. 


6. Calculation of V's, V/¥YLWL's and normalizing V/YLWL 

me Vos and ۷/1 و‎ are calculated based on the Vain >» Vincrement 
LWL, and m calculated above. The ٦/۴59 are then normalized using the 
uation in Table II. If Vrax = Vmin 2nd Vincrement = 0., only one speed 


is used. 


IE. Calculation of m for a 400 foot ship and PC for a 600 foot ship 
The C, for a ship of L= 400 ft. (Cep) and PC for a ship of L = 600 
ft. are calculated for each V/YLWL- using the normalized parameters determined 


| x 
above, and equations (7) and (8) respectively. Subroutine CEF is used to 





۲ See program lisving in Appendix A. 
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evaluate the Cnebyshev Polynomials for these equations. 


€. Calculation of C for a 400 foot ship and Cp for a ship of length L 


The Reynolds Numbers (Re) for the L= 400 ft. ship (Re ) and for the 
ship of length L (Re.) are calculated for each V/YLWL by using equations 
(11) and (12). 


pe” - x v. m) 
۱۳ ۱ x 1052 


: (12) 
Re u LWL x 3 


5 
v 
The C.'s for the hOO ft. ship (Cp, ) and ship of length L (Cra) are 


calculated using the ATTC 1947 correlation line given in equation (13). 


cage 


f = 


1.132 log, (Re x Cp) (13) 


In the calculation of the Ce رگ‎ it is necessary for the program to use an 
| iteration process, since 7 appears on both sides of equation (13). The 
| ۰۰۰٠ eszimate Tor this process is made by using the equation for the 


| International Towing Tank Conference (ITIC) line, equation (14). 


ch 9-97 > (14) 
3E 3 
981g Re 2) 
Equation (13) is recycled for 50 times, or until the Co values converge to 
Within a tolerance of 5.0 x 10 EC. Although this should not occur, و‎ 


formula (13) does not converge, the program prints the error signal 


mote that LWL = 406.75 tor a ship of 400 ft. length between perpendiculers, 


Ca 
and that kinemetie viscosity is for salt water at 59 5. 


ES 





|"SCHOENHER NO CONVERGZ" with the speed at which the trouble was experienced. 


9. Calculation of wetted surface (S) if not input 

really, the wetted surface (S), in square Teer ISTE +٣ Lane 
computer as input. However, if the wetted surface is not known at that 
Game, the program will calculate an approximate wetted surface using 


equation (15). 


۱ -.06 . 666 | 
S = 2,654 (T/B) pes (B/H) ES (©) L (V) (15) 


This equation for approximating wetted surface for the Series 60 was 
developed by D. Khoushy in reference (2). This formula was tested E 
Sively and compared to the contours for wetted surface calculation, Figures 
METAL through B-126, in reference (1). It was founi Wes thie formula 
adequately represents m wetted surface contcurs of the Series 60, with 
the maximum error rarely exceeding 1.2%. Since wetted surface is a factor 
ecu tines resistance, EHP and SHP, if wetted surface is not an input 
CT co the program, one can expect errors in EHP and SHP due to the 
use of equation (15), in addition to the errors to be expected due to the 


use of equations (7) and (8) as discussed in Chapter III. 


mo. Calculation of = ee Ship of length L, R, EHP, and SHP 
The e for the ship of length L, (Cr), is calculated using equation 
6 (۰ 


pa ے + ںاد‎ + ۷٥ (16) 
ts tr fu fs f 


Then the remaining calculations are trivial, namely: 


ES 








EP OE S ٤۶ 





ts 
EHP = R v 
S 
550 
SHP = EHP 
ES 


Mere PC is for a ship of L= 600 ft. 


EN Calculation of LCB position from as a percentage of L 
The LCB position from 3 as a percentage of length between perpen- 

ers is calculated and presented in the output. It should be —— 
EE CESt corrections to R, EHP, ahd SHP for variation in LCB وه‎ from 
those of the parent model hull 0 are not built into the program. Graphs 
mor assistance in making these corrections, Vasea Ou eS 49 through 53 

of reference (1), ere included in Appendix B. The LCB position from yf 

EN percentage of lengtn between perpendiculars is calculated for entry 


Lo these graphs. 


l2. tout 
The output for each ship calculated includes وت‎ L/B, B/H, Cp? Cm, 
length between perpendiculars (LBP) in feet, LWL in feet, B in feet, Hin 


feet, displacement (DISP) in tons, volume of displacement (VOL) in cubic 


feet, wetted surface (WSURF) in square feet, Alo, water density P in 


2 


ID. secó , kinematic viscosity (GNU) in ft^ , LCB position as a percentage 


— | sec 

of LBP, and a tabulation of V, V/YLWL, R, EHP, PC, and SHP for each speed 
requested in the input. A statement is printed under this tabulation which 
Meninds the user trat the values are not corrected for variation in LCB 


position as explained above. 


دو 





۷ و‎ EXAMPLES IN USE OF THE EHP-SAP COMPUTER FROGRAM FER KNOWN SHIPS 


fee timed EXHP-SHP computer program was tested with the parameters 
ES Six ships whose EHP-SHP characteristics are known. As a first test, 
the Series 60 009 of the ship SCHUYLER OTIS BLAND was chosen, 
since this was the ship used in Appendix D of reference (1) in illustrating 
a numerical example of the use of the Series 60 charts. It was also 
Sen since it illustrates the ability of the program to interpolate for 
intermdeiate values of = (0.651) and L/B (6.82) and to extrapolate 
outside the range of the Series 60 in B/H (2.414). 

As a second test, the five parent models of the Series 60, with | 


L = 600 feet, were chosen since complete information on their EHP and SHP 


D 


behavior is zresent2d in Tables 27 through 31 in reference (1), and they 
cover a variety of hull proportion combinations. 


Inner results of these tests are discussed below. 


۱ test using the ship SCHUYLER OTIS BLAND 

As a first test on this ship, input Mode I was used with the wetted 
surface area of 39994.0 square feet as input. As a result, the computer 
output, Table VI, for EHP was a measure of the accuracy of equation (7). 
mumconparing the computer program's EHP results with those resulting from 
the hand calculations carried out in Table D-3 of reference (1), it was 
found that the maximum error experienced in EHP was + 2.70% at ۷ھ‎ = .85, 
while the average error experienced was + 1.19% . Although the average error 
of 1.19% is fairly gratifying, the maximum error certainly is not. However, 
referring to Table IV again we see that in using the hull proportions of the 


CHUYLZR OTIS BLAND we place ourselves in the neighborhood of parent model 
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4218 where the maximum error of 2.74% was experienced in the test of 
equation (7). 

Tne SCHUYLER OTIS BLAND was also tested without wetted surface as 
MG. ine results of this run are presented in Table VII. We can see 
that the wetted surface calculated by equation (15), 39651.54 square feet, 
1s in error by -.858%. The error in the calculation of wetted surface area 
Bends to reduce the error in EHP in this instance. 

One must keep in mind that we &re comparing computer output with 
Wealues derived from hand calculations, in which design charts were read with 
an unknown accuracy. 

2. Test using the five parent models of Series 60 with L = 600 feet 

estilos of the computer run with the tive parent models of 
Series 60 with L = 600 feet are presented in Tables Vill through XII. In 
these cases, tne wetted surface area of the ships was calculated by equation 
B5), since it سی‎ known. The maximum and average percent errors in 
ERP EN c crprenesdeun this test are presented in Table XIII. Although 
at some points the errors exceeded 2%, in no case did the average error 
exceed 1%. These errors are probably due mainly to the innacuracy of the 


equation used for wetted surface area. 





VI. CONCLUSIONS AND RECOMMENDATIONS 


ESESg upon the results of tests run to date, one can conclude that 
ee EHP-SEP computer program presented herein can be relied upon to 
estirate the EHP-SHP requirements of single-screw merchant ships, falling 
within the range of hull proportions of the Series 60, to a tolerance in the 
neighborhood of 1% on the average, and with great certainty, to within a 
tolerance of less than 3% for any given speed. 
The weak link in the EHP-SHP program appears to be the calculation 
x 때 0۵4 surface area. It is recommended that the representation of the 
Series 60 wetted surface contours be further investigated, to determine 
whether equation (15) can be modified and improved. 
it is furtner recommended that the elements of matrices [a] and 
ES 11 be punched in the form of FORTRAN data statements, and included in the 
body of the main program. This will greatly facilitate its conversion to 
ENcubroucine for inclusion in a larger program if it were so desired. 
meio elso recommended that the program be extensively tested with 


EUusting ships, to further determine its accuracy. 
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APPENDIX A 


PROGRAM LISTINGS 


۱ Flow chart for EHP-SHP program. 
2. Listing of EHP-SHP program and S/R CER. 


۲ listing of the elements of matrices JA and | 8 š 
nl nl 
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۱1۲۱ 17 37-7 
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ar 
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Input deck for matrices A ajend fB] 


Mie input of ratrices fas} and ERE a total of 900 elements, 
is accomplished by reading them from a deck of 180 cards, with 5 
elements on each card, arrengsd in Format 5F16.7. See Appendix A for 
their proper arrangenent. 


Total no. of cards = 180 


eee identification card 

In columns 1 - 4, (NJOBS), Punch the integer number of ships to be 
calculated in this particular run, with the 
units digit of the integer in column k. 
Format Ih. 


fin columns 5 - 80, (QQ) Punch any identification information desired. 
Example: "SCHUYLER OTIS BLAND TEST RUN." 
Format 19A. 


Total no. of cards = 1 


ee Input for Snips 

a. If input Mode I is being used for a ship, punch the following 3 cards 
mor each ship: 

Card 1 


In columns 1 - 4, (INOPT), Punch a positive integer. Example: "+5". 
Format Il. 


Card 2 


In columns 1 - 10, (XLBP), Punch the length between perpendiculars in feet. 
Example: "400.0" Format F10.2. 


In columns 11 - 20, (BM), Punch the maximum beam at the load water line in 
feet. Example: "66.0". Format F10.2. 


In columns 21 - 30, (H), Punch the draft in feet. Example: "27.0". 
Format: F10.2. 


In columns 31 - HO, (CB), Punch the block coefficient. Example: "0.651". 
Roura ۲۳1022 | | 


In columns 41 - 50, (VMIN), Punch the lowest speed in knots for which you 


ama LION OA Sue and SAP. Example: 
ST. Format F10.2 


-'[5- 





In columns 51 - 60, (VMAX), Punch the highest speed in knots for which 
penwdesire iniorration on EHP and SHP. 
Example: "19.25". Format F10.2. 


In columns 61 - 70, (VINC), Punch the increment of speed at which you wish 
option On EHP and SHP. Example: "1.07" 
OEE FIO.2. 


Card 3 


In columns 1 - 10, (DCF), Punch the frictional ship correlation factor 
you desire to use. Example: ".oook". 
Forrat F10.2. 


In columns ll - 20, (XLCB), Punch the position of LCB from WÉ in feet, 
with a+ sign if forward and a - sign if 
aft of 8 . Examoles “- 420". Format F10.2. 


In columns 21 - 30, (WSURF), Punch the wetted surface area in square feet. 
Example: "39994.0". Format F10.2. If wetted 
surface is not known, punch a "0.0" in this 
space and the program will calculate it. 


In colurns 5. - +0, (P=) Punch the density of water in lb. sec“ /£t. 
at the temperature desired. Example: "1.9905". 
Format PEOS Z: 


In columns 41 - 50, (GNU), Punch the kinematic viscosity in ft.“/sec. 
Example "1.281-05", for 1.2817 x 107. 
Porras ElO.2. | 


Ec ImhPorsxbcs0Hn OM ENP and SHP is desired for only one speed for the 
TD نهد‎ 1۳۳ speed, on card 2 punch this speed in the spaces for 
VMIN and VMAX, and punch "0.0" in the space for VINC. 





Toll No. OI cards = NJOBS X 32 
۳ lupe mode Ti is being used for a ship, follo: the instructions 
690. 0460 © above mach the following exceptions: 
Card Ll 


In columns 1 - 4, (INOPT), Punch a negative integer. Example: "-5". 
Format Ih. 


Card 2 


In columns 11 - 20, (BH), Punch the beam to draft ratio. Example : 
۱ Maso Format F10.2 


In columns 21 - 30, (DEL), - Punch the displacement in tons. Example: 
"10000.0". Format F10.2. 


E 





In columns 31 - LO, (CP), Punch the prismatic coefficient. Example 
mar, Format F10.2. 


Card 3 
figenticel to Mode I. 
Mode I and Mode II input ships can be mixed at will as long as the 
cards are punched as above. 
Graphs for making corrections CO CMe computer output tor variation 


CB position ere contained in Appendix B. 


سد 


4. Loading of cards 
Ine Cars Oa an tie following order: 
J. Mein EHP-SHP rrosram Deck 
2-۰۰۲517567017 61016 CEF 
5.۰7060 0" ۲0 matrices la.) anc [3,1] 
h. The identification card 


3-۳۷۲۰ ۲٦۲36 ۲٦0" 0ت6‎ ۰ 


= = 





